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We suggest and demonstrate a tomographic method to fully characterize homodyne detectors at the quantum 
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validate the reconstructed operator measure on nonclassical states. Finally, we exploit results to estimate the 
overall quantum efficiency of the detector. 
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I. INTRODUCTION 

Balanced homodyne detection is a crucial detection tech¬ 
nique for continuous variable quantum technology and lies at 
the core of many quantum optics experiments I EE). From its 
first proposal for the measurement of quadrature squeezing, 
to its current extensive use in the fields of quantum tomogra¬ 
phy, quantum communication and quantum metrology gm, 
this detection scheme has carved its place into experimental 
quantum optics. Besides quantum optical systems, homodyne 
detection extends its reach to the whole field of continuous 
variable quantum technologies, spanning from atomic systems 
|[T6l 17]] to quantum optomechanics tm. 

Advances in technology promoted the spread of many dif¬ 
ferent configurations of this versatile apparatus, tailored to 
disparate experimental needs. Such a wide range of appli¬ 
cations calls for reliable ways to fully characterize homo¬ 
dyne detectors. Each specific setup relies on classical cali¬ 
brations in order to gain the most general description of the 
apparatus and of the relationship between the input state and 
the measurement output. Recently, a characterisation of ho¬ 
modyne detection, only used as a phase-insensitive photon 
counter, was demonstrated CD. However, a general and re¬ 
liable model for the description of the fully phase-sensitive 
homodyne detection in the form of a Quantum Detector To¬ 
mography (QDT) is, in fact, still lacking. 

The pioneering proposals for QDT (201423) were followed 
by the experimental characterization of an avalanche photo¬ 
diode, in both single and time-multiplexed configurations, for 
the detection of up to eight photons (24l . Subsequent works 
developed the idea, including the effect of decoherence onto 
the operator description (25), or different detection devices, 
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such as superconducting nanowires (26) and TES based sys¬ 
tems (221 2 Jl 1 . However, these are detectors devoted to pho¬ 
ton counting, whose description is entirely embedded in the 
diagonal sector of the Fock space. Only recently, a specific 
phase-sensitive hybrid scheme, in the form of a weak homo¬ 
dyne detector based on photon counting, was the object of 
an experimentally realized QDT [ 29 ]. In this paper, we move 
several crucial steps forward, and present a theoretical and ex¬ 
perimental realization of QDT for homodyne detector, i.e. the 
most commonly used form of a fully phase-sensitive detector, 
whose operators are naturally described in phase space. 

The quantum description of any detector is given by a pos¬ 
itive operator-valued measure (POVM), i.e. a set of positive 
operators {n n }, giving a resolution of identity n n = I. 
The determination of these operators is, in turn, the main goal 
of detector tomography. Given an input state p, the Born 
rule states that p £ = Tr[pn n ] is the probability of obtain¬ 
ing the outcome n when the generalized observable described 
by n is being measured. The inversion of this formula allows 
the reconstruction of the operators n n from the experimen¬ 
tally sampled probability distribution p p n , over a suitable set of 
known states p. These must form a tomographically complete 
set, spanning the Hilbert subspace where the POVM elements 
are defined on Eo). 

The simplest choice for a continuous variable system is a 
set of coherent states. They provide an overcomplete basis 
for the Fock space, and it has already been proved that even 
1-dimensional discrete collections of coherent states form a 
complete basis, and may be used to reconstruct classical and 
non-classical states (3TJ 32]. In fact, the experimental distri¬ 
butions of the outcomes for a set of coherent states already 
provide a full representation of the detector operators, in the 
form of a sample of their Q-functions 

Q n (a) = -Hll n |a) = -P“, 

7T 7T 

where {Vf} represents the probability distribution for a co¬ 
herent state. However, this representation is not suitable to 
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provide a complete and reliable characterization of the detec¬ 
tor. In fact, any subsequent use of this reconstruction scheme 
to predict the outcome of the measurement for a different sig¬ 
nal would involve the (numerical) evaluation of the trace rule 
in the phase-space as 

Tr[pII n ] = [d 2 aP[p](a)Q n (a), 

Jc 

where the Glauber P-function P[p\ ( a ) is singular for any non- 
classical state, and thus not suitable for sampling. In order 
to overcome this problem we suggest an expansion in the 
quadrature basis of the operator measure associated with the 
detector, using as probe a set of coherent states. We then ob¬ 
tain the coefficients of the expansion using a least squares al¬ 
gorithm on a sufficiently large sample of data. We also val¬ 
idate the experimentally obtained POVM by reconstructing 
nonclassical known states. Finally, we exploit results to es¬ 
timate the overall quantum efficiency of the detector. 

The paper is structured as follows. In Section [IT] we re¬ 
view the description of homodyne detection and introduce the 
algorithm employed for the reconstruction of its POVM. In 
Section [lll| we describe our experimental apparatus, whereas 
in Section|IV|we present results of the reconstruction, as well 
as their validation on nonclassical states. Section [V] closes the 
paper with some concluding remarks. 

II. HOMODYNE DETECTION 

A homodyne detector is a fully phase sensitive apparatus 
that provides a complete characterization of any given state 
of a single-mode radiation field l33l . This state, the signal , 
is sent to a balanced beam splitter, where it interferes with 
an intense coherent field, the local oscillator , usually coming 
from the same laser source. The phase of the signal has then a 
precise value </> with respect to the local oscillator, and can be 
adjusted by means of a piezo-actuated mirror. The two outputs 
of the beam splitter are then focused on two photodiodes, and 
the resulting photocurrents subtracted and analyzed. It can be 
shown that, in the approximation of high amplitude \/3\ of the 
local oscillator, the measurement associated to this detector 
corresponds to 

atf+a^b + 

—2— a) 

where a and b are the mode operators for the signal and the 
local oscillator, respectively. The operator b was replaced by 
\/3\ in Eq. {!]) by considering its action onto the local oscil¬ 
lator, that can be treated as a coherent state \/3). The operation 
connected to the working scheme of this detector is therefore 
the measurement of the quadrature operator x$ on the signal 
mode. Such a link states the equivalence between the dis¬ 
crete spectrum of the operator X and the continuous one of 
the quadrature, due to the high intensity of the local oscillator, 
that can be consequently treated classically. This equivalence 
can be extended to the characteristic functions 

Tr a , b [e lXx p 0 |/3) </?|] ^ Tr a [e lXx *p] (2) 


assuring the equivalence of all moments El. 

The key feature of homodyne detection is its ability to dis¬ 
cern between different phase values of an input signal, set¬ 
ting itself apart from the photon counting detectors that have 
been characterized in the past. A straightforward choice for a 
basis in which representing the POVM elements of a phase 
insensitive device is the number basis, in the form 11^ = 
\ n )( n \* where \n)(n\ is the projector onto the n- 
photon Fock state. Such a description is no longer suitable 
for our apparatus, that hinges on a phase-sensitive operation 
scheme. Off-diagonal elements in the number-basis expan¬ 
sion could enclose phase-sensitive properties, as was done in 
|[29lL but reconstruction of the detector operators would be¬ 
come increasingly difficult due to the high dimension of the 
Hilbert space the POVM would be defined in. 


A. The reconstruction algorithm 

A suitable basis to expand the POVM H(x) of a phase- 
sensitive detector for continuous variable systems is the set 
{\y)(y\} y eR °f eigenstates of a quadrature operator e.g. 
setting (j) = 0 

U(x) = Jdyg { y x) \y){y\ . (3) 

The detector has now a Q-function representation given by 

Va{x) = (a|n(x) |a) 

= J d V9 ( y x) \J^ exp {— 2 (y ~ M) 2 } > (4) 

(x) 

where the matrix g k ’ contains all the information needed to 
describe the response of the detector. The set of operators 
{n(x)} can then be discretized, reflecting the experimental 
sampling during the measurement process, and be confined to 
a selected portion of the quadrature range, say [x m i n , x max ]. 
The expansion on the quadrature basis {\y)(y\} can be dis¬ 
cretized as well, reducing the number of POVM elements. 
Equation ([4]) may be rewritten as 

n i = n (^) = Yl 9 k \vk){vk\ • (5) 

k 

The analogue of Eq. ([4]), i.e. V a (xj) = (a\U(xj)\a), can 
now be compared to experimental results and may be inverted 
to find the matrix g. To this purpose, starting from a set of 
coherent states with calibrated amplitudes {|a s |}, we use a 
least-squares method: 

= arg min i ^ ^ \p &s (xj) - V &s (xj)] 2 l , (6) 

H}>o { 3 s J 

where Pa s (%j) is the experimentally observed distribution for 
the coherent state with amplitude \a s \. A pictorial represen¬ 
tation of the algorithm is presented in Fig. [TJa)- The algo¬ 
rithm retrieves the matrix g by comparing the experimentally 
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sampled distributions to the corresponding Q - function rep¬ 
resentation. Starting from a coherent state of amplitude \a s \ 
and a quadrature value Xj , the algorithm finds the set of coef¬ 
ficients g k that minimises the difference between (xj) and 
Va s (xj). Experimental imperfections and fluctuations, cor¬ 
responding to the shaded blue area in Fig. [TJa ), degrade this 
relationship and “switch on” new coefficients in the expansion 
of Eq. 0- At the same time, each calibrated amplitude \a s \ 
is associated to a rescaled value |a s | = |a s | j s , to explicitly 
include the detector finite efficiency. The least squares algo¬ 
rithm of Eq. 0 performs this minimisation simultaneously 
for all the quadrature values Xj G [x m { n , x max ] and all the 
coherent states in the set. The detector response is then fully 
contained in the matrix g. 



Figure 1. (Color online) Quantum tomography of a homodyne de¬ 
tector. A pictorial representation of the least-squares algorithm is 
reported in panel (a). The experimentally sampled probability distri¬ 
butions pa ( x ) (solid blue) are compared to the Q - function repre¬ 
sentation V&(x) (dashed yellow). This is shown dashed to visualise 
the discretisation of Eq. (3]). Panel (b) shows the efficiency of re¬ 
construction as a function of increasing amplitude spacing in the to¬ 
mographic set. The function A (5a) is here normalized to the matrix 
dimension N . A sharp transition can be seen as the spacing between 
the states of the set approaches the maximum for the overlap func¬ 
tion of Eq. 0- Panel (c) shows the phase space representation of the 
full set of coherent states used for tomography. Panel (d) shows the 
result of the reconstruction obtained imposing | a s \ = \a s \. Notice 
the (expected) diagonal shape. 

If we look back at Eq. 0, it is quite natural to link the char¬ 
acteristics of the matrix g to the features of the detector recon¬ 
struction. Each matrix row g k associates a quadrature value Xj 
to a set of projectors on quadrature eigenstates, with weights 
given by the coefficients in Eq. 0. For an ideal detector, the 
matrix is diagonal g k = 5k j , i.e. the only nonzero coefficient 
associates a quadrature value Xj to its projector \xj)(xj\. As 
we have seen, experimental imperfections will degrade this 
one-to-one relationship, spreading the coefficients around a 
central value. At the same time, the full array of coefficients 
7 S provides a unified model for the response of the detector. A 
detector tomography devised in this way is general enough for 


application to different configurations of the homodyne detec¬ 
tion. 

Let us now focus on the tomographic set, and in particu¬ 
lar on the characteristics required to perform a reliable recon¬ 
struction. To this aim we have performed simulated experi¬ 
ments with sets having an increasing number of equidistant 
coherent states, with amplitudes in a given range a G [—3, 3], 
and measured at the same phase 0 = 0. For each set the matrix 
representation of the detector g k was retrieved, and associated 
to the amplitude spacing 5a between the coherent states in the 
set. Since the identity matrix I is the ideal-case solution of the 
reconstruction algorithm, we consider the following function 
of the amplitude spacing 5a 


A (5a) = Tr I 



(7) 


as a figure of merit to assess tomographic sets, e.g. to find the 
minimal 5a corresponding to a reliable reconstruction. 

For a perfect reconstruction, the two matrices are both the 
identity and A reduces to the dimension of the matrix g. In the 
opposite case, more and more elements on the diagonal will 
be voided, and the trace will decrease. In Fig. |TJb) we show 
the results: a steep transition, corresponding to a deterioration 
of the reconstruction, appears for 5a ~ 0.7. A similar conclu¬ 
sion may be obtained theoretically upon considering the over¬ 
lap of two Gaussian p a (x, xo) of equal standard deviation a 
but varying center xq, i.e. 


/ +oo 

dxp a (x,xo)pa(x,x 0 + y ). 

-OO 


( 8 ) 


In particular, the function y • f(y) may be used to assess the 
tomographic set of coherent states, as it captures, roughly 
speaking, the trade-off between an increasing spacing and a 
decreasing overlap. Upon substituting a = 0.5, as it is for 
coherent states, we have that y • / (y) has a maximum at 
I/max = 2 -1 / 2 0.707, in good agreement with the value 

obtained by simulated experiments via Eq. 0. The reliabil¬ 
ity of this estimation has been then confirmed experimentally 
(see below). 


III. EXPERIMENTAL APPARATUS 

A schematic diagram of the experimental setup is shown 
in Fig. [2] The apparatus is based on a mode-locked 
Ti:sapphire laser (Spectra-Physics Tsunami) providing, after 
suitable splitting, both the local oscillator (LO) beam for bal¬ 
anced homodyne detection and the probe coherent states for 
detector tomography. The laser emits 2 ps pulses at a cen¬ 
tral wavelength of 785 nm, with a repetition rate of 82 MHz. 
The detector characterized in this paper is an optical homo¬ 
dyne apparatus, operating in the time domain at high sampling 
frequency f35j|36]|. Amplitudes of the probe coherent states 
were selected by means of a reflective-coating glass attenu¬ 
ator. Precise calibration of each state is done by means of 
a Type I BBO crystal cut for degenerate spontaneous down- 
conversion (SPDC), pumped by the frequency doubled por¬ 
tion of the main laser beam. The injection of the probe coher¬ 
ent states into the signal path of the SPDC triggers the stim¬ 
ulated emission of downconverted photon pairs in the same 
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mode (thus generating single-photon-added coherent states, 
SPACSs 1373) and in the idler mode, at a rate proportional 
to 1 + | a | 2 , where a is the amplitude of the incoming coher¬ 
ent state. Such a procedure provides a precise standard-free 
calibration of the input amplitude lf38l by means of the ratio 
between the count rate of stimulated and spontaneous events 
pursuing the idea of a calibration-free characterization. 



Figure 2. (Color online) Schematic diagram of the experimental 
setup. Picosecond-duration pulses at 785 nm are split by two high- 
transmission beam splitters. A first beam serves as the local oscilla¬ 
tor (LO) for balanced homodyne detection (B.H.D.), while a second 
provides the probe coherent states for detector tomography. Opti¬ 
cal delay lines provide synchronization, while the probe amplitude 
is adjusted with the attenuator (A). The main portion of the beam is 
frequency doubled in a LB O crystal, and then spatially and spectrally 
filtered by means of a pin-hole and a Michelson interferometer (F). 
It finally serves as the pump for a type-I BBO crystal for paramet¬ 
ric down conversion, whose idler mode is used for calibration means 
and, after narrow Etalon spectral filtering (E), as a trigger signal. 


IV. DETECTOR TOMOGRAPHY 


In our detector tomography we have focused attention to 
the range x E [—2, 2] probed by a tomographic set of coher¬ 
ent states with amplitudes a s E {—3, 3}. Our set is made of 
coherent states with 12 different amplitudes, each one mea¬ 
sured at 9 different phase values between 0 and n. The full set 
is represented in Fig.[ljc). 

As a first step we have performed a preliminary validating 
step of our reconstructing algorithm, by neglecting the rescal¬ 
ing factors 7 S . The set of amplitudes {a s } has been measured 
with the homodyne detector, as the mean value of the proba¬ 
bility distributions, and the coefficient matrix g has been re¬ 
trieved imposing \a s \ = \a 3 \. The result, reported in Fig.[ljd), 
presents the expected diagonal shape, partially blurred by fluc¬ 
tuations. The homodyne detection model obtained in this way 
may efficiently describe the detector behavior on several clas¬ 
sical and non classical states. 

Upon confirmation of the effectiveness of our method, we 
selected a reduced collection of coherent states, concentrat¬ 
ing on those with a null or i r phase difference between the 


signal and the local oscillator. This set of 25 coherent states 
has been then inserted into the algorithm of Eq. ([ 6 ]) together 
with the set of calibrated amplitudes { a s }, and the results are 
shown in Fig.[3ja). The retrieved matrix showed the expected 
diagonal shape, with an experimental spread around the cen¬ 
tral value. The array of coefficients is reported in Fig.[3jb). 
An estimate for the uncertainty is obtained from the inverse of 
the second derivative of the function minimized in the least- 
squares algorithm. This uncertainty is increasing for decreas¬ 
ing values of a, since for the limiting case of a = 0 our model 
is not defined. It was then used to weight the , and get a final 
mean value of 70 = 0.90 =b 0.03. 



Figure 3. (Color online) Quantum tomography of a homodyne de¬ 
tector. Panel (a) shows the matrix g J k representing the expansion of 
the POVM elements in the quadrature basis. By comparison of this 
matrix with that of Fig. 0 d) one notices the spread of the coefficients 
around the diagonal, due to the additional parameters 7 S . These are 
reported in blue in the plot of Fig. [3jb), where the shaded area rep¬ 
resents an estimate of the uncertainty, obtained through the second 
derivative of the function in Eq. In a similar fashion, the mean 
value 70 and its uncertainty are indicated in solid and shaded yellow. 
These results were then validated on the reconstruction of nonclas- 
sical states. In panels (c) and (d) we report the experimentally re¬ 
constructed Wigner functions for the Fock state |1) and for a single 
photon-added coherent state. The side plots in both panels show the 
Wigner marginal distributions (yellow), which are both very close to 
their representations in the detector description (the uncertainty due 
to the spread of 70 is given by the blue-shadowed areas). 


A. POVM validation 

The aim of quantum detector tomography is to fully char¬ 
acterize a given detector by retrieving the set of operators 
that fully describes its measurement process. In order to 
validate our technique for QDT, we have employed the re¬ 
constructed POVM to reproduce measurements performed on 
known states. In particular, we have tested the POVM on two 
kinds of nonclassical states: a single-photon Fock state and a 
SPACS with a = 0.5. These two states have been first mea- 
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sured in our setup, and the experimental results have been then 
compared to those obtained for the same states using the re¬ 
constructed detector POVM. The photon addition scheme of 
our setup, fundamental for the generation of these two states, 
had been previously characterised (39), and the operator a t 
had been found to apply to a given input state with preparation 
efficiency ( « 0.91. The measured quadrature distributions 
for the Fock state and the SPACS, reported in Fig. [3jc) and 
[3jd), are in excellent agreement with the expected Q-function 
representation based on the tomography of our homodyne de¬ 
tector. The robustness and reliability of our method has been 
thus confirmed and we proved that the specific experimental 
realization of the detector, which depends on several parame¬ 
ters (like the detrector quantum efficiency, the degree of mode 
matching, the alignment, etc.), can be efficiently captured by 
the tomographic procedure. We also proved that the results of 
subsequent measurements can be effectively reproduced. 



Figure 4. (Color online) Changing the number of coherent states in 
the tomographic set. Panels (a) and (b) show results of a success¬ 
ful tomography for the limiting case of an amplitude spacing of 0.7, 
as predicted by the theoretical modeling. The limiting case of over- 
sampling was also tested, yielding another positive result. Results 
of QDT using the full tomographic set of 109 states are reported in 
panels (c) and (d). 


The tomographic set of coherent states that we have used 
throughout the paper has proved to be a good test for our de¬ 
tectors. In order to improve the accuracy of the reconstruc¬ 
tion we may extend the set to cover a bigger portion of phase 
space, while to minimize the experimental effort we may want 
to reduce the number of states in the set. Proceeding as above, 
and following our predictions from Eq. ([7]), we found that 
a minimal set of nine coherent states may be selected from 
the experimental data, all at phase zero. The amplitude spac¬ 
ing is three times larger than in the previous situation, but the 
reduced set is still able to provide a quorum for the tomog¬ 
raphy. Results are presented in Fig. Qa) and (b). On the 


other hand, even with the full set of 109 coherent states we 
have been able to efficiently reconstruct the detector, despite 
the increased phase space coverage and the increased fluctua¬ 
tions, due to oversampling. The matrix g° k and the coefficients 
7 s for this case are shown in Fig. |4|c) and (d). 

On the basis of the previous analysis and considering the 
large uncertainty of data in the area close to the origin, we 
have assigned a fixed value = 1 to probe states with 
amplitudes smaller than 0.5. We found that the values of 7 
for the smallest and the largest tomographic set are given by 
71 = 0.91 =b 0.03 and 72 = 0.84 zb 0.09 respectively. In 
fact, all the values of 70, 71 and 72 are comparable within the 
uncertainty. On the other hand, they convey different informa¬ 
tion regarding the detector. The increased coefficient spread 
of the matrices reported in Fig. [3ja) and Fig. |4|a) can be 
considered as an additional rescaling parameter, modeling the 
experimental fluctuations, that is therefore directly included in 
the tomography. Values for 7 are then larger, with reduced un¬ 
certainty. The matrix of Fig. |4jc) has instead a smaller spread, 
and therefore the extra rescaling is conveyed in 72 , lowering 
its value and making it more accurate, even though less pre¬ 
cise. Thus 72 is better suited to be compared to the value of 
the quantum detection efficiency that can be obtained by clas¬ 
sical calibration. Indeed, we find that 72 ~ y/fj, despite the 
fact that our model does not involve any prior knowledge of 
the detector structure or implementation. Different sets can 
therefore be used to highlight specific properties of the detec¬ 
tor, adding value to our technique. 


V. CONCLUSIONS 


We have suggested and demonstrated a full quantum de¬ 
tector tomography technique for a homodyne detector. In 
ideal conditions each detector operator is associated to a single 
quadrature projector: our technique suitably describes how 
experimental noise and specific physical realizations of the 
detector affect this description and allows us to quantify ex¬ 
perimentally the spreading of the detector operators onto ad¬ 
jacent quadrature states. The model is general enough to de¬ 
scribe any kind of homodyne setup, and it has proven capa¬ 
ble of effectively describing the detector response to different 
tomographic sets. The reconstructed POVM have been then 
validated on different nonclassical states, thus confirming the 
robustness and the reliability of the method. 

Our results provide a general method to estimate the overall 
detection efficiency in this class of detectors and may repre¬ 
sent a valuable resource to optimize homodyne detection in 
different situations. Our model may be generalized to specifi¬ 
cally treat single parameters of homodyne detectors, as mode 
mismatch or correlations between amplitude and phase noise. 
Besides, a better understanding of the fundamental function¬ 
ing of this detector paves the way to an evolution of the same, 
as well as a broader and more precise use in quantum optics 
and quantum technology with continuous variables. 
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